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Abstract 

We study the inclusive hadrodroduction of D^, D^, D*^, and Df mesons at next-to-leading 
order in the parton model of quantum chromo dynamics endowed with universal non-perturbative 
fragmentation functions (FFs) fitted to e~^e~ annihilation data from CERN LEPl. Working in the 
general-mass variable-flavor-number scheme, we resum the large logarithms through the evolution 
of the FFs and, at the same time, retain the full dependence on the charm-quark mass with- 
out additional theoretical assumptions. In this way, the cross section distributions in transverse 
momentum recently measured by the CDF Collaboration in run II at the Fermilab Tevatron are 
described within errors. 
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Recently, there has been much interest in the study of charmed-hadron (Xc) production at 
hadron colhders, both experimentally and theoretically. The CDF Collaboration measured 
the differential cross sections da/dpT for the inclusive production of D^, D^, D*^, and Z)+ 
mesons (and their antiparticles) in pp collisions at the Fermilab Tevatron (run II) as functions 
of transverse momentum {px) in the central rapidity (?/) region jl| . Unfortunately, the most 
advanced theoretical predictions available so far 3, Isj, based on quantum chromodynamics 
(QCD) at next-to-leading order (NLO), consistently undershoot all the D^, D^, and D*^ 
data by significant amounts, as is evident from Fig. 3 of Ref. Q], while no predictions for 
mesons exist yet. It is presently an open question if this discrepancy is related to an 
experimental problem, a technically deficient QCD prediction, or the appearance of new 
physics beyond the standard model. Such a situation is familiar from inclusive bottom- 
flavored-hadron (Xt) production in run I, where a long-standing discrepancy between CDF 
data and certain NLO predictions of QCD were, in fact, interpreted as an indication for 
low-energy supersymmetry [7]. It is, therefore, an urgent task to deepen our understanding 
of the inclusive hadroproduction of charmed hadrons on the basis of QCD in order to render 
the theoretical predictions as reliable as possible, so as to establish a sturdy anchor for 
new-physics searches. This is even more important in view of future physics at the CERN 
Large Hadron Collider, where the continuum production of charmed hadrons will provide 
important backgrounds for numerous new-physics signals. This is the main motivation of 
this letter. Moreover, we provide the first NLO prediction for the CDF data Prior 
to explaining our improved theoretical framework and describing our updated input, for the 
reader's quick orientation, we present a brief survey of the various NLO approaches adopted 
so far in the literature. 

In the so-called massless scheme, also known as zero-mass variable-flavor-number (ZM- 
VFN) scheme 0,0], the conventional parton model approach implemented in the modified 
minimal subtraction (MS) scheme is adopted, assuming that the charm (c) quark can be 
treated as massless, although its mass m is certainly larger than the asymptotic scale pa- 
rameter Aqcd- In this approach, the c quark occurs as an incoming parton, leading to 
contributions in addition to those where it is produced by an incoming gluon g or a. light 
u, d, or s quark. The c quark fragments into the charmed hadron similarly as the gluon 
and the light quarks with a fragmentation function (FF), which must be known from other 
processes. The well-known factorization theorem provides a straightforward procedure for 
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systematic higher-order perturbative calculations. Due to the assumption that m = 0, the 
predictions are reliable only for large values of px, with px ^ m, where powers of m? /p^ 
can be neglected. This approach has the advantage, however, that the potentially large log- 
arithms of the type ln(p|,/m^) are absorbed into the c-quark parton distribution functions 
(PDFs) of the colliding hadrons and into the FF for the transition c X^. These logarithms 
are resummed through the Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) evolution 
equations. Still, m appears in the initial conditions of the c-quark PDF and FF. In this 
respect, the c-quark PDF and FF differ from the PDFs of the gluon and the light quarks. 

In the so-called massive scheme, also called fixed-flavor- number (FFN) scheme Q], the 
number of active quark flavors in the initial state is fixed to = 3, and the c quark only 
appears in the final state. The physical value of m is explicitly taken into account together 
with the variable pt, as if the two were of the same order. In this scheme, m acts as a cutoff 
for the initial- and final-state collinear singularities. However, in NLO, terms proportional to 
\n{p'^/rn?) arise at large pt values from collinear gluon emission off c quarks or from almost 
collinear branchings of gluons into cc pairs. For px ^ m, these terms spoil the convergence 
of the perturbative series. The FFN approach with = 3 is thus limited to a narrow pt 
range, reaching up to a few times m. The advantage of this scheme is that the m^/p^ power 
terms are fully taken into account. 

Obviously, the ZM-VFN and FFN schemes are valid in complementary p^ regions, and is 
desirable to combine them in a unified approach that enjoys the virtues of both schemes and, 
at the same time, is bare of their flaws, i.e. one that resums the large logarithms, retains 
the full finite-m effects, and preserves the universality of the FFs. This is vital for a reliable 
and meaningful interpretation of the CDF data ^1, which mostly populate the transition 
region. A first attempt to implement such an interpolating scheme is the so-called fixed- 
order next-to-leading-logarithm (FONNL) scheme, in which the conventional cross section 
in the FFN approach is linearly combined with a suitably modified cross section in the ZM- 
VFN approach with perturbative FFs, using a suitable pr-dependent weight function Q . 
In both finite-mass approaches, FFN and FONNL, the FFN cross sections are convoluted 
with a non-perturbative c-quark FF, adjusted to e^e^ data, that is not subject to DGLAP 
evolution. 

Here, we wish to advocate an approach that is much closer in spirit to the ZM-VFN 
scheme, but keeps all m^/p|. power terms in the hard-scattering cross sections, namely the 
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general-mass variable-flavor-number (GM-VFN) scheme, which has recently been elaborated 
for the photo- and hadroproduction jl^ . Q] of heavy-flavored hadrons. In this approach, 
one starts from the ^ m region and absorbs the large logarithms ln(p|./m^) into the c- 
quark PDF of the incoming hadrons and the FF for the c — > Xc transition. After factorizing 
the Inm^ terms, the cross section is infrared safe in the limit m — 0, and = 4 is taken 
in the strong-coupling constant ag and the DGLAP evolution equations. The remaining m 
dependence, i.e. the m? /p^ power terms, is retained in the hard-scattering cross sections. 
These terms are important in the intermediate pt region, where pt ^ fn, and are expected 
to improve the precision of the theoretical predictions. The large logarithms are absorbed 
into the PDFs and FFs by subtraction of the collinearly (mass) singular terms. However, 
in order to deflne a unique factorization prescription, one also has to specify non-singular 
terms. This is done by requiring that, in the limit px oo, the known ZM-VFN hard- 
scattering cross sections are recovered. To achieve this, subtraction terms are derived by 
comparing the FFN theory in the limit m — with the ZM-VFN theory, implemented in 
the MS factorization scheme. This matching procedure is useful, since all commonly used 
c-quark PDFs and FFs are deflned in the ZM-VFN scheme. The latter can then be used 
consistently together with hard-scattering cross sections calculated in the GM-VFN scheme. 
The derivation of the subtraction terms is described in Ref. [isj ]. 

We note that our implementation of the GM-VFN scheme is similar to the Aivazis- 
Collins-Olness-Tung (ACOT) jl^ scheme formulated for the initial state of fully inclusive 
deep-inelastic scattering. The extension of this scheme to the inclusive production of heavy 
partons was considered in Ref. where the resummation of the flnal-state coUinear log- 
arithms was only performed to leading order (LO) and parton-to-hadron FFs were not 
included. A comprehensive discussion of the differences between our approach and Ref. ^| 
concerning the coUinear subtraction terms may be found in Ref. jisl . 

We now describe our calculation of the differential cross section dP a / [dpTdy) of p + p — *• 
Xc + X, where Xc = D^, D^, D*^, and X comprises the residual flnal-state hadrons, at 
NLO in the GM-VFN scheme. A crucial ingredient entering this calculation are the non- 
perturbative FFs for the transitions a X^ where a = g,u, u, d, d, s, s, c, c. For Xc = D*^, 
such FFs were extracted at LO and NLO in the MS factorization scheme with Uf = 5 
massless quark flavors several years ago ^ from the scaled-energy {x) distributions da/dx 
of the cross section of e+ + e" ^ D*+ + X measured by the ALEPH and OPAL 17 1 
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Collaborations at CERN LEPl. Two of us recently extended the ^ analysis of Ref. ^ to 
include = D^,D^,D^,A^ by exploiting appropriate OPAL data [l9]. Besides the total 
Xc yield, which receives contributions from Z ^ cc and Z ^ bb decays as well as from light- 
quark and gluon fragmentation, the ALEPH and OPAL Collaborations separately specified 
the contribution due to tagged Z bb events yielding X^ hadrons, which then weakly decay 
to Xc hadrons. The contribution due to the fragmentation of primary c quarks into Xc 
hadrons approximately corresponds to the difference of these two measured distributions. 
To test the scaling violations of these FFs and also the separation of the c ^ X^ component, 
these FFs were employed to interpret the x distributions of e"^ + e~ center-of- 
mass energy 75 = 10.55 GeV measured by the CLEO Collaboration J2l| at LEPP CESR, 
with very encouraging results. Further details may be found in Ref. |l 

In Refs. 0, the starting scales /iq for the DGLAP evolution of the a — »• Xc 
FFs in the factorization scale taken to be /io = 2m, with m = 1.5 GeV, for 

a = g,u,u,d,d,s,s,c,c and /io = 2mb, with nib = 5 GeV, for a = b,b. The FFs for 
a = g,u,u,d,d,s,s are assumed to be zero at fi'p = and are generated through the 
DGLAP evolution to larger values of fi'p. Since the effect of the gluon FF is important at 
Tevatron energies, as was found for D*^ production in Ref. we repeated the fits of the 
Xc FFs for the choice fio = ttl, nib. This changes the c-quark FFs only marginally, but has 
a strong effect on the gluon FF. For shortage of space, these new FFs will be presented 
elsewhere. 

The calculation of the cross section d'^a/{dpTdy) proceeds as outlined in Ref. The 
full cross section consists of three contributions. The first one contains all the channels 
with only gluons or light quarks in the initial state and c-quark fragmentation. Only this 
contribution carries explicit m dependence. Second, this contribution must be extended by 
allowing for c quarks in the initial state. The third contribution is due to gluon or light- 
quark fragmentation. The second and third contributions are calculated in the ZM-VFN 
scheme using the hard-scattering cross sections derived in Ref. [2^. A certain part of these 
contributions is due to Feynman diagrams with internal c-quark lines; another one is due 
to diagrams with external c-quark lines and contains m-dependent logarithms, which are 
resummed. Thus, in the FFN scheme, the m dependence of these contributions would only 
enter beyond NLO, which is refiected in the ZM-VFN scheme by the generic suppression of 
the c-quark PDF relative to the gluon and light-quark ones and of the gluon and light-quark 
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FFs relative to the c-quark one. This entitles us to omit this m dependence by calculating 
the c-quark-initiated contributions and those involving the fragmentation of gluons or light 
quarks in the ZM-VFN scheme. It turns out that the light-quark fragmentation contributions 
are negligible. However, gluon fragmentation contributes approximately 40% to the cross 
section, almost independent of px- For the D*'^ case, we showed in Ref. ^l^] that the effect 
of the m-dependent terms is much reduced in the full cross section, since those parts that 
have to be calculated with m = dominate. In fact, this observation carries over to the 
other charmed mesons considered here. 

We are now in a position to present our numerical results for the cross sections of inclusive 
D^, D^, D*^, and -D^ hadroproduction to be directly compared with the CDF data Q|, 
which come as distributions da/dpT at = 1.96 TeV with y integrated over the range 
\y\ < 1. For each species, the particle and antiparticle contributions are averaged. We 
work in the GM-VFN scheme with nj = 4, thus excluding hadrons from Xf,-hadron 
decays, which are vetoed in the CDF analysis We set m = 1.5 GeV and evaluate 



22l | , corresponding to 



tts"'^'*(/i_R), where fj.R is the renormalization scale, with = 328 M{ 
a^^\mz) = 0.1181. We employ proton PDF set CTEQ6.1M from the CTEQ Collaboration 



22| | and the FFs introduced above. We distinguish between the initial- and final-state 



factorization scales, fip and fi'p, so that we have three unphysical mass scales altogether. 



Our default choice is fiR = fip = f^'p = where rriT = y^T + transverse mass. 

In order to conservatively estimate the theoretical error due to the scale uncertainty, we 
independently vary the values of fin/rriT, ^p/fnT, and fip/rriT between 1/2 and 2, and 
determine the maximum upward and downward deviations from our default predictions. 

Our theoretical predictions are compared with the CDF data on an absolute scale in Fig.^ 
and in the data-over-theory representation with respect to our default results in Fig. |21 The 
four frames in each figure refer to D'^, D~^, D*~^, and mesons. In all cases, we find 
good agreement in the sense that the theoretical and experimental errors overlap, i.e. the 
notorious discrepancy between experiment and theory jl| mentioned in the introduction has 
disappeared. In fact, our theoretical predictions provide the best description of the CDF 
data obtained so far. 

As for the D'^, D*~^, and mesons, many of the central data points fall into the 
theoretical error band, while those for the mesons lie somewhat above it. With the 
exception of the case, the experimental results are gathered on the upper side of the 



6 




FIG. 1: Comparison of the CDF data ]l] with our NLO predictions for = Z)°, D+,D*+,D+. The 
sohd lines represent our default predictions, while the dashed lines indicate the scale uncertainty 
(see text). 
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FIG. 2: Data-over-theory representation of Fig. ^with respect to our default predictions. 



theoretical error band, corresponding to a small value of fiR and large values of fip and 
fi'p, the ^Ft dependence being dominant in the upper pt range. As is evident from Fig. |21 
in these cases, the central data points tend to overshoot the central QCD predictions by a 
factor of about 1.5 at the lower end of the considered pt range, where the errors are largest, 
however. This factor is rapidly approaching unity as the value of pt is increased. The 
tendency of measurements of inclusive hadroproduction in Tevatron run II to prefer smaller 
renormalization scales is familiar from single jets, which actually favor /ir = pt/2 ^2^. 

The overwhelming bulk of the theoretical error stems from the scale uncertainty discussed 
above. Residual sources of theoretical uncertainty include the variations in the value of m 
and the adopted PDF and FF sets. We now quantitatively study the impact of these 
variations relative to the typical example of our default prediction for D*~^ mesons. The 
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generous variation of m by ±20%, from 1.2 to 1.8 GeV, induces a shift in cross section of 
only ±2% at = 5 GeV, which rapidly decreases towards larger values of pt because the 
m-dependent terms are themselves reduced in size. Switching to the NLO proton PDF set 
MRST2004 of Martin, Roberts, Stirling, and Thorne 



2^, with = 347 MeV and a more 



physical parameterization of the gluon distribution to produce a better description of the 
Tevatron inclusive jet data, produces a reduction ranging from —15% at = 5 GeV to 



-1% at = 20 GeV. The theoretical uncertainty due to the FFs was estimated in Ref. 



to be of order 10% or less in the pt range considered here, 
separately to slightly incompatible ALEPH and OPAL 



w comparing FF sets fitted 
1711 data. 
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In conclusion, the GM-VFN scheme, which we elaborated at NLO for the inclusive photo- 
[lll | and hadroproduction jl^ . \\\ of heavy-flavored hadrons, resums large logarithms by the 
DGLAP evolution of non-perturbative FFs and guarantees the universality of the latter as in 
the ZM-VFN scheme and simultaneously retains the m-dependent terms of the FFN scheme 
without additional theoretical assumptions. Adopting this framework in combination with 



Id 



3, 



we managed for the 



new fits of D+, D*+ , and FFs to OPAL data from LEPl 
irst time to reconcile the CDF data on the production of these mesons in Tevatron run II 
3] with QCD within errors and thus eliminated a worrisome discrepancy. Furthermore, we 
presented the first NLO predictions for the data [l|. 
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